Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



Intermediate resolution H/3 spectroscopy and photometric 

monitoring of 3C 390.3 

I. Further evidence of a nuclear accretion disk * 

' A.I. Shapovalova^, A.N. Burenkov^, L. Carrasco'^''^, V.H. Chavushyan-^ , V.T. Doroshenko'*, A.M. Dumont^, 
■ V.M. Lyuty^ J.R. Valdes^, V.V. Vlasuyk\ N.G. Bochkarev^ S. Collin^, F. Legrand^, V.P. Mikhailov^, 

O.I. Spiridonova^, O. Kurtanidze^, and M.G. Nikolashvili® 



o 

(N 



^ Special Astrophysical Observatory of the Russian AS, Nizhnij Arkhyz, Karachaevo-Cherkesia, 369167, Russia 
e-mail: ashapOsao.ru 

. Instituto Nacional de Astrofi'sica, Optica y Electronica, Apartado Postal 51, CP. 72000, Puebla, Pue., Mexico 

e-mail: carrasco@inaoep .mx 
^ Observatorio Astroriomico Nacional, UNAM, Apartado Postal 877, CP. 22860, Ensenada B.C., Mexico 
* Sternberg Astronomical Institute, University of Moscow, Universitetskij Prospect 13, Moscow 119899, Rusia 



o 



> 

' ^ DAEC, Observatorie de Paris, section de Meudon, Place Jansen, F-92195 Meudon, France 

, ^ Abastumani Astrophysical Observatory, Georgian AS, Mt. Kanobili, 383762, Abastumani, Georgia 

\D 

O ! Received: February 21, 2001 / Accepted: June 13, 2001 

O I Abstract. We have monitored the AGN 3C 390.3 between 1995 and 2000. A historical B-Band light curve dating 

, back to 1966, shows a large increase in brightness during 1970 - 1971, followed by a gradual decrease down to a 

^ l' minimum in 1982. During the 1995-2000 lapse the broad H/3 emission and the continuum flux varied by a factor of 

« 3. Two large amplitude outbursts, of different duration, in continuum and H/3 light were observed ie.: in October 
1994 a brighter flare that lasted «1000 days and in July 1997 another one that lasted a; 700 days were detected. 
^ The response time lag of the emission lines relative to flux changes of the continuum, has been found to vary with 

, time ie. during 1995-1997 a lag of about 100 days is evident, while during 1998-1999 a double valued lag of ~ 100 

days and ~ 35 days is present in our data. The flux in the H/3 wings and line core vary simultaneously, a behavior 
indicative of predominantly circular motions in the BLR. Important changes of the H/3 emission profiles were 
, detected: at times, we found profiles with prominent asymmetric wings, as those normaly seen in Syls, while at 

» I ' other times, we observe profiles with weak almost symmetrical wings, similar to those of Syl.8s. We further dismiss 

' the hypothesis that the double peaked H/3 profiles in this object originate in a massive binary BH. Instead, we 

found that the radial velocity difference between the red and blue bumps is anticorrelated with the light curves 
of H/3 and continuum radiation. This implies that the zone that contributes with most of the energy to the 
emitted line, changes in radius within the disk. The velocity difference increases, corresponding to smaller radii, 
as the continuum flux decreases. When the continuum flux increases the hump velocity difference decreases. These 
transient phenomena are expected to result from the variable accretion rate close to the central source. The optical 
continuum and the H/3 flux variations might be related to changes in X-ray emission modulated by a variable 
accre tion rate, changing the surface temperature of the disk, as a result of a variable X-ray irradiation ( Ulrich, 



2000). Theoretical H/3 proflles were computed for an accretion disk, the observed proflles are best reproduced by 
an inclined disk (25°) whose region of maximum emission is located roughly at 200 Rg. The mass of the black 
hole in 3C 390.3, estima ted fr om the reverberation analysis is Mrev ~ 2.1 x 10^ Mq, 5 times larger than previous 



estimates (Wandel et al. 199!: ) 



Key words, galaxies: active - galaxies: Seyfert - galaxies: individual (3C 390.3) - line: profiles 



Send offprint requests to: A.I. Shapovalova or L. Carrasco 

* based on observations carried out at SAO RAS, Guillermo 
Haro, Abastumani and Sternberg Observatories. 



1. Introduction 

Quasars and Active Galactic Nuclei (AGNs) are amongst 
the most luminous objects in the Universe and provide 
a unique test case for several processes of astrophysical 
importance. A relevant characteristic of some of these ob- 
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jects is a variable broad emission line spectrum. The region 
where broad lines are formed (hereinafter BLR) is rather 
close to the central engine and may hold basic informa- 
tion to the understanding of the formation and fueling of 
AGNs and quasars, a process not yet fully understood. 

Long temporal baseline spectral monitoring of the nu- 
clei of some AGNs has revealed a time lag in the response 
of the broad emission lines relative to flux changes of 
the continuum. This lag depends on the size, geometry 
and physical conditions of the emitting region. Thus, the 
search for correlations between nuclear continuum changes 
and flux variations of the broad emission lines may serve 
as a tool for mapping the geomet rical and dynamical 
structure of the BLR (see Peterson, 1993 and references 
therein). From the study of the responses of lines of dif- 
ferent ions, to changes of the nuclear continuum and the 
comparison of these data with photoionization model pre- 
dictions, one could infer some of the physical conditions 
in the BLR as a functi on of the distance from the central 
source (Dumont et al. 1998). 

From changes of the broad emission line profiles, we 
could distinguish, in principle, different kinematical mod- 
els for the BLRs, including cases in which matter is under- 
going accelerated motions, ie. falling or outflowing, keple- 
rian rotation in the gravitational held of the central body. 



and some others (Bochkarev & Antokhin 1982; Blanford 
& McKee |1982| Antokhin & Bochkarev |1983|). 



During the last decade, the study of the BLRs in some 
objects has met considerable success, mainly due to an 
increasing number of coordinated multiwavelength moni- 
toring campaigns through the international "AGN Watch" 
program. This program has provided data with good tem- 
poral frequency and coverage for a num ber o f selected 
Seyfe rt galaxies : NGC 5548 (Clavel et al.|l99l| ; Peterson 
et al.|l99l| , |l993| , |l994 |l999|; M aoz et al. |l993i Die trich et 
al. | 199^ ; Korista et al.|l995| ; Chiang et al. ^OOO|) ; NGC 



et al. 
al. 



3783 (Reic hert et al. |l994 Stirpe et a l |1994| ; AUoin 



1995|) ; NGC 4151 (Cre nshaw et al. |l996|; K aspi et 



199^ ; Warwick et al. |l996t Ed elson et al. |l996| ); Faira ll 



9 (Rodriguez-Pascual et al. 1997 ; Santos-Lle o et a l. 1997) 
NGC7469 (Wanders et al |l997t Colher et al.|l99^ ) and 3C 
390.3 (Leighly et al. |l997|; Dietrich et al. |l998|; O'Brien et 



al. 1998). It has been found that a significant part of the 



BLR variability can be associated with a region located 
within some light weeks from the central source. The in- 
fered BLR extent is comparable to the size of an hypothet- 
ical accretion disk surrounding a supermassive black hole. 
Therefore, at least part of the flux from the BLR arises, 
apparently, from the accretion disk itself (Laor & Netzer 
1989|, Dumont & Collin-Souffrin |19904 |1990b| ; Zheng et 



al. 1991; Hubeny et al. 2000). It is also found that, the 



higher ionization lines respond faster to continuum flux 
changes than those of lower ionization and that the op- 
tical and ultraviolet continua vary without a significant 
time lag between them. 

The results of the studies of the velocity fields in the 
central sources are still ambiguous. It is not clearly estab- 
lished yet, whether or not, the time lags in the response 



of the blue and red wings and the core of the line pro- 
file with respect to continuum variations are the same. So 
far, studies of the velocity dependent response of the C iv 
A1549 emission line in NGC 5548 have found no evidence 
of significant radial motions in its BLR. (Korista et al. 



1995) 



Most of the objects, included in the multiwavelength 
monitoring AGN Watch, are radio-quiet Syl galaxies and 
only one, 3C 390.3 (z=0.0561), is a well known broad line 
radiogalaxy. This type of object represents about 10% of 
the radio-loud AGNs. It is also a powerful double lobed 
FRII radio-galaxy with a relatively strong compact core. 
The two extended lobes, at a position angle of 144°, each 
one with a hot spot, are separated by about 223"(Leahy 
& Perley 1991). A faint well-collimated thin jet at P. A. = 
37°, connecting the core to the no rthern lobe, has been ob- 
served by Leahy & Perley (|l995| ). The VLBI observations 
at 5 GHz show evidence of super luminal moti on (w ith v/c 
--4) in this parsec scale jet (Alef et al. |1996| , [19881) . 

The strong variability of this object, both in the con- 
tinuum light and the emission lines is well known (Barr 
et al. 119801; Yee & Oke |l98l|; Netzer |l982l Barr et al. 



1983; Pension & Perez 



Veilleux & Zheng 1991 



1996; Wamsteker et al. 



1997; Dietrich et al. 1998; O'Brien 



1984 Clavel & Wa mstek er |l 987 
Shapovalova et al. 1996| ; Zheng 



et al. 1995). The object is also a highly variable X-ray 



source, with a spectrum showing a broad Fe Ka line (Inda 
et al. |1994|; Eracleous et al. |l996t Wozniak et al. |1998D. 



During a multiwavelength monitoring campaign in 1995, 
several large-amplitude X-ray flares were observed in 3C 
390.3; in one of them the X-ray flux increased by a factor 



of 3 in 12 days (Leighly et al. 1997). Evenmore, Leighly 
& O'Brien (1997) have presented evidence for nonlinear 
X-ray variability. 

In an analysis of lUE spectra of 3C 390.3, obt ained 
during the 1978-1986 period, Clavel & Wamsteker |l987 



detected a variability time lag of 50 and 60 days between 
the broad emission lines of Civ A1549 and Lya, and the 
UV-continuum. However, from a rean alysis of the same 
spectral data, Wamsteker et al. 1997 have derived a lag 



of 116 ± 60 days for Civ and 143 ± 60 days for Lya. 
Furthermore, from lUE monitoring data for December 
1994 to March 1996, O'Brien et al. (|l998| ) obtained a cor- 
responding lag of 35 days for the Civ A1549 emission line 
and 60 days for Lya. Yet, from optical monitoring of 3C 
390.3 in the 1994-1995 period, Dietrich et al. ( |1998[ ) de- 
rived a time lag of about 20 days for the Balmer lines re- 
sponse to changes in the X-ray continuum. Furthermore, 
at that time, no temporal lag between the optical and 
the UV or X-ray continua changes was detected. The UV- 
bump, usually observed in a large number of Seyfert galax- 
ies is very weak or even absent in 3C 390.3 (Wamsteker 
et al. 1997). Our object is a prototype of a class that 
shows very broad, double peaked emission line profiles, to 
explain this profile, a number of possible scenarios have 
been suggested. These models contemplate the following 
hypothetical physical models: 
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Supermassive binary black holes (Gaskel 1983). If each 
of the black holes (BH) had its own associated BLR, 
then the blue and red bumps, seen in the broad lines, 
could be assesed to the corresponding BH. As a re- 
sult of precession of the binary BH, secular displace- 
ments of the peak wavelengths are expected. From the 
analysis of 3C 390.3 spectra, obtained during 20 years 
(1968-1988), Gaskel ( |1996D detected a displacement of 
the blue bump in H/3 at -1-1500 km/s, and determined 
from the radial velocity changes a possible period of 
about 300 years that corresponds to 7x 10^ Mq for the 



1998 



1 999 



2000 



BHs. However, Eracleous et al. ( 1997 ), having added 
data on the displacement of the Ha blue peak for the 
1988-1996 lapse, obtained a period for the binary BH, 
P > 800 years, with an estimated total mass for the 
binary BH larger than 10^^ Mq. These authors argued 
that such a large value for the mass of the BH is diffi- 
cult to reconcile with the fact that, the BHs found so 
far in the centers of galaxies, have masses lower than 
10^ Mq. Therefore, they reject the hypothesis of a bi- 
nary black hole for 3C 390.3. 

Outflowing biconical gas streams (Zheng et al. 1991). 
In this case the response of the broad emission lines to 
variability of the continuum should occur from the blue 
wing to the red wing passing through the line core, ie. 
the blue and red wings must show different time lags. 
Although Wamsteker et al. (1997) claimed that the 
Lya and Civ A 1549 blue wings lagged behind the red 
ones, implying radial gas motions, this result is not 
confirmed by Zheng ( 19961 ). Apparently, these contra- 
dictory results are caused by the rather poor temporal 
sampling of the UV data (~94 days). Dietrich et al. 
( 1998D have not detected any lag in the response of the 
blue and red wings relative to the core of Ha and H/? 
emission. It was also found that the flux in the wings 
of H/3, Lya and C i v A15 49 vary quasi-si multa neously 
(Shapovalova et al. |l996| ; O'Brien et al. |l998|) . These 
results imply that there are no significant radial mo- 
tions in the BLR, and Livio & Xu (1997) have shown 
that the double-peaked lines seen in 3C 390.3 cannot 
be produced by two line-emitting streams, since the 
emission region of the far cone jet would be obscured 
by the optically thick accretion disk. 
Broad Line emission from an accretion disk in 3C 
390.3 has been suggested by Perez et al. (1988). In 
this case, the emission from an inclined (30°) accre- 
tion disk provides a good fit to the broad Ha and 
H/3 profiles. (Rokaki et al. 1992). Simple disk models 
predict two symmetrical bumps, the blue bump be- 
ing slightly brighter due to doppler-boosting effects. 
However, lUE data shows a blue wing of Lya much 
brighter than the red one, this cannot be explained by 
a doppler boosting effects alone. There are also epochs 
when the red wing of Lya is brighter than the blue 
one (Zheng 1996), i.e. the flux ratio of the red to the 
blue wings is quite variable. Also, quasiperiodical vari- 
ations {P « lOyr ) of the observed flux ratios of the 
blue and red wings of H/3 have been reported (Veilleux 
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Fig. 1. The fight curves of 3C 390.3 in the B,V and Rc 
photometric bands for a 10 "aperture, in the 1997-2000 
period. 



1991; Shapovalova et al. 199(]| ; Bochkarev et 



& Zheng 

al. 1997a). In order to explain these variable ratios, 
more complex models are required; for instance, hot 
spots (Zheng et al. 1991); two-arm sp iral w aves in the 
accretion disk (Chakrabarti & Wiita 1994 ); or a rela- 
tivistic eccentric disk (Eracleous et al. 1995). 
Long temporal base line studies of the changes in the 
broad emission line profile through systematic moni- 
toring shall allow us to distinguish between those mod- 
els. 

In this paper, we present the results of the spectral 
(H/3) and photometric BVRI monitoring of 3C 390.3 
during the 1995-1999 period. This work is part of the 
long-term monitoring program for about 10 Seyfert 
galaxies of different nuclear luminosities started in 
1986 at the SAO RAS, and carried out jointly, since 
1998, with the INAOE (Mexico), and several ob- 
servatories in the Former Soviet Union (FSU) and 
West-European countries (Bochkarev et al 
Bochkarev & Shapovalova 1999| ). 



1997b 



2. Observations and Data Reduction 
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2.1. Photometry 



Table 2. Photometric Results 



The BVR photometry of 3C 390.3 was carried out at 
three different observatories i.e. SAO RAS in North 
Caucasus (Russia), the Crimean Laboratory of the 
Sternberg Astronomical Institute (CL SAI, Russia) and 
the Abastumani Astrophysical Observatory (AAO) in 
Georgia. The observations at the SAO RAS during 1997- 
2000 were obtained with the 1 m and 60 cm Zeiss tele- 
scopes equipped with an offset guided automatic photome- 
ter. The instrument has a liquid nitrogen cooled CCD 
camera with a format of 1040 x 1160 pixels, (Amirkhanian 
et al. 2000). The scale at the CCD being 0.45"per pixel, 
with a corresponding field of view of 7.5 x 8.5 arcmin. 
Exposures of the morning and evening sky were adopted 
as flat-field frames. Bias and dark current frames were also 
obtained. Data reduction was carried out with the soft- 
ware developed at SAO by Vlasyuk ( |1993| ). The photome- 
try is obtained by signal integration in concentric circular 
apertures, of increasing size, centered at the baricentcr of 
the measured object. The instrumental photometric sys- 
tem of this instrument is similar to those of Johnson in B 



Table 2 is only available in electronic form at the CDS via 
anonym ous ftp to cdsarc.u-strasbg.fr (130.79.125. 5) or via 
jittp:/ /cdsweb. u-strasbg.fr/ Abstract. htmi 



and V, and of Cousins (1976) in R and I. The VRI pho- 



tometry during 1998-1999 was obtained at the Crimean 
Laboratory of the Sternberg Astronomical Institute with 
the 60 cm telescope (Zeiss-600) equipped with an ST-6 
thermoelectrically cooled CCD camera. The photometric 
system of CL SAI is similar to that of Johnson. Further de- 
tails about SAO an d Crim ean observations can be found in 
Doroshenko et al. (|200l|) . BVRI observations at the AAO 
in Georgia during 1997-1998 were obtained with the 70 cm 
menisk telescope equipped with a CCD camera. These 
data were reduced with the IRAF package (DAOPHOT). 
As local photometric standard, we adopted stars "B" and 
"D" of Pension et al. (|l97l|) which are close to 3C 390.3 in 
our CCD images, consequently, effects due to differential 
air mass are negligible for internal calibration purposes. 
The adopted BVRI magnitudes for the reference stars are 
given in Table |l|. 

Table 1. Adopted Photometry for the Standard Stars 



Star 


B" 


ya 




I*" 


B 


15.04 


14.28 


14.13 


13.59 


D 


15.40 


14.65 


14.42 


13.90 




"From Penston et al. 
From Dietrich et al. 



The results of the broad-band photometry of 3C 390.3 in 
the BVRc filters for a circular aperture of 10"obtained 
at AAO and CL SAI were transformed to the SAO BVR 
photometric system. Our results are presented in graph- 
ical form in Figurej^. The estimated mean errors in the 
BVRI photometry for the entire data set are 0.026, 0.021, 
0.020, 0.022 magnitudes, respectively. In Table || our pho- 
tometric results and the corresponding errors are listed. 




4000 4500 5000 5500 

Wavelength [A] 

Fig. 3. The spectra corresponding to the high activity 
(top) and the low activity (bottom) states 

The last column lists, the observatory where the data was 
obtained i.e. (SAO, SAI — Crimean Laboratory, Abast. - 
Abastumani Astrophysical Observatory in Georgia) . 

2.2. Optical Spectroscopy 
2.2.1. Observations 

Spectra of 3C 390.3 were obtained with the 6 m and 
1 m telescopes of the SAO RAS (Russia, 1995-1999) 
and at INAOE's 2.1 m telescope of the "Guillermo Haro 
Observatory" (GHO) at Cananea, Sonora, Mexico (1998- 
1999). The spectra were obtained with long slit spectro- 
graphs, equipped with CCD detector arrays. The typi- 
cal wavelength interval covered is that from 4000 A to 
7500 A. The spectral resolution varied between 4 and 15 
A. Spectrophotometric standard stars were observed ev- 
ery night. The specific information about the instrumental 
set-ups in different telescopes is listed in Table ||, there we 
list: 1 - telescope; 2 - type of focus; 3 - spectrograph; 
4 - CCD format; 5 - set-up code. The log of spectro- 
scopic observations is given in Table I listed in Cols. 1 - 
UT date; 2 - Julian date; 3 - code according to Table 3; 
4 - projected spectrograph entrance apertures (the first 
dimension being the slit-width, and the second one the 
slit-length); 5 - wavelength range covered; 6 - spectral 
resolution; 7 - seeing; 8 - signal to noise ratio in a region 
of the continuum (5370-5420A) , where there are no promi- 
nent emission or absorption lines. The spectrophotometric 
data reduction was carried out either with software devel- 
oped at SAO RAS by Vlasyuk(1993) or with the IRAF 
package for the spectra obtained in Mexico. The image 
reduction process included bias, flat-field corrections, cos- 
mic ray removal, 2D wavelength linearization, background 
subtraction, stacking of the spectra for every set up, and 
flux calibration. 
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Fig. 2. A historical light curve for 3C 
390.3. Open cicle s - Sandage (1973), 
Neizvestny (1986); large filled circles 
- photometry (present paper ); open 
squares - Yee & Oke (1981); small 
filled circles - Ba badzhanyan ts et al. 



1984) 



crosses 
Scott et al. 
(|197^), Pica et 



( [1973| , |1974 [1975| , |1976 
- Can non et al. ( 1971 ^ 
( [1976|), S elmes et al 
al. (1980); plusses - spectral contin- 
uum (present paper); open triangles - 
Dietrich et al. (199q); filled triangles- 
spectral cont. from Perez et al. (1988) 
and Lawrence et al. (1996). 



Table 3. Characteristics of the telescopes and spectrographs 



Telescope 


Focus 


Equipm. 


CCD 
(pixels) 


Code 


1 


2 


3 


4 


5 


6m SAO 


Prime 


UAGS 


530x580 


G 


6m SAO 


Prime 


UAGS 


1024x1024 


P 


6m SAO 


Prime 


MPFS^ 


530x580 


T 


6m SAO 


Nasmyth 


Long slit 


1024x1024 


N 


Im Zeiss 


Cassegrain 


UAGS 


530x580 


Zl 


Im Zeiss 


Cassegrain 


UAGS 


1040x1170 


Z2 


2.1m GHO 


Cassegrain 


B&Ch 


1024x1024 


M 



MPFS refers to a Multi-Pupil Field Spectrograph 



Table 4. Log of the Spectroscopic Observations 

Table 4 is only available in electronic form at the CDS via 
anonym ous ftp to cdsarc.u-strasbg.fr (130.79.125.5 ) or via 



http: / /cdsweb.u-strasbg.fr/Abstract.html 



2.2.2. Absolute Calibration of the Spectra 

Since, even under good photometric conditions, the accu- 
racy of spectrophotometric measurements is rarely better 
than 10%, for the study of AGN's variability the standard 
scheme for flux calibration by means of comparison with 
stars of known spectral energy distribution is not accept- 
able. Instead, we use the flux of the narrow forbidden emis- 
sion lines which in AGNs are non- variable in timescales of 
tens of years. This is largely due to the extension of the 
narrow line emitting region (NLR) and to the lower gas 
density present in it. Thus, the effects of light traveling 
through this region, jointly with the long recombination 
times (troc)~100 years, for Uo ^^10"^ cm~^), damp out the 



short time scale variability. Hence, the bright narrow emis- 
sion lines are usually adopted as internal calibration for 
scaling AGN spectra (Peterson 1993). However, in the case 
of 3C 390.3 there has been some discussion about the vari- 
ability of these lines. Let us address this question in more 
detail. 



From lUE data, Clavel & Wamsteker (|1987D found that 
the flux in the narrow components of Lya and C IV A1549 
emission lines in 3C 390.3 decreased continuously by « 
40% between 1978 and 1986. The involved timescales sug- 
gest an upper limit of about 10 light years for the size of 
the regions in which these lines are formed. 



Zheng et al. ( 1995| ) reported a considerable variabil- 
ity of the [Oiii] A4959 narrow emission line fluxes (~1.8 
times) from observations carried out between 1974 and 
1990. We believe that this result is somewhat doubtful 
in the light of the following considerations: the variability 
is only infcrcd when comparing results from observations 
carried out with different instrumental set-ups. They in- 
clude, on the one hand, fluxes from image dissector scan- 
ner (IDS) data, obtained with the Shane 3 m telescope at 
Lick Observatory through an aperture of 2.8 x 4" (1974- 
1983 lapse) and on the other hand, observations carried 
out at La Palma with the 2.5 m Isaac Newton telescope 
(INT), equipped with an image photon counting system 
(IPCS) through a 1.6" wide slit (1984-1989). The low- 
est flux values largely correspond to those observed at La 
Palma through this rather narrow slit. Hence, it is pos- 
sible that the large differences in flux reported for the 
[O III] A4959 line could be due to effects associated with 
the narrow slit used. Evenmore, in the spectra published 
by Zheng et al. ( 1995 ), it is clearly seen that, during 1988 
and 1989, the red wing of H/3 showed an important en- 
hancement and its contribution could have also affected 
the determination of the [O ill] line flux. Unfortunately, 
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these authors do not provide a detailed discussion of their 
adopted procedure for flux determination. 

Relevant in the context of forbidden line variability are 
the observations by Yee & Oke (1981) that covered the 
1969-1980 lapse, these were carried out with the Hale 5 m 
telescope using a multichannel spectrophotometer through 
an aperture of ~10". They found that the [Oiii] lines 
and the narrow Balmer line fluxes did not changed during 
those years. Furthermore, a detailed study of the [O ill] 
A5007 flux variability in 3C 390.3 was carried out by 
Dietrich et al. ( |1998D during 1994-1995. The spectra were 
obtained through wide apertures (5"x7.5"and 4"xl0") 
and were calibrated by comparison with broad band pho- 
tometry. It was found, as expected, that the [O ill] A5007 
line flux remained constant (within 2.6%) during the ob- 
served time interval. Therefore, we consider that in the 
case of 3C 390.3 there is no reliable evidence of the [O ill] 
AA4959-I-5007 flux variability in time scales of months to 
few years. Our spectra for the H/3 region were scaled by 
the [O III] AA4959-I-5007 integrated hue flux under the as- 
sumption that the latter did not change during the time 
interval covered by our observations (1995-1999). A value 
of 1.7x10"" ergs s'^cm'^ (Veilleux & Zheng |199lD for 
the integrated [O III] line flux was adopted. In order to 
calculate a normalization coefficient, the continuum was 
determined in two 30 A wide clean -line free- windows 
centered at 4800 A and 5420 Arespectively. After contin- 
uum subtraction, blend separation of the H/3 and [O ill] 
components was carried out by means of a Gaussian fit- 
ting procedure, applied to the following: H/3 — broad 
blue, broad red and central narrow; [Olll] AA4959,5007 
— broad and narrow components. The forbidden lines are 
represented by two Gaussian curves with an intensity ratio 
I(5007)/I(4959)=2.96. 



2.2.3. Intercalibration of the Spectral Data 

The flux of H/3 was determined from scaled spectra after 
continuum subtraction and removal of the [Oili] doublet. 
The narrow component of H/3 was not subtracted due to 
possible ambiguities in the fitting scheme for the broad 
and narrow components of the line. We should mention, 
that the contribution of the narrow component to the in- 
tegral H/3 flux is only about 7%. As Peterson & Collins 
( |1983| ) pointed out, it is necessary to correct the fluxes 
for aperture effects, due to the fact that while the BLR 
and nonstellar continuum are effectively pointlike sources 
(<1"), the NLR is an extended one (>2-15"). Then, the 
measured NLR flux depends on the size of the spectro- 
graph's entrance aperture, and in the case of direct images, 
on the size of the adopted aperture. The light contribution 
of the host galaxy to the continuum depends also on the 
aperture size. Peterson et al. ( 1995D have shown that the 
uncertainties in CCD photometry can be minimized, when 
sufficiently large apertures are adopted. For an aperture 
of 5"x7.5", the expected photometric errors are typically 
2 to 3 %. 



The NLR in 3C 390.3 is more compact than in most Syl 
galaxies, in narrow-band [O ill] images, this object shows 
a compact nuclear emission without signs of an extended 
structure (Baum et al. 19881 ). The results of panoramic 



two-dimensional spectrophotometry of the nuclear region 
of this object show that the [Oiii] A5007 emission arises 
from a zone smaller than r < 2"(Bochkarev et al. 1997a). 
Furthermore, Osterbrock et al. ( |1975| ) obtained a rather 
low value for the [Oiii] F(A4363)/F(A5007) line ratio, im- 
plying a moderately high electron density in the NLR (sev- 
eral lO^cm^"^). Since the NLR in 3C 390.3 can be consid- 
ered as a point source, we did not applied corrections for 
aperture effects to the nonstellar continuum to narrow line 
flux ratios or to the broad to narrow line flux ratios, as 
the light losses in the slit are similar for these compo- 
nents. However, the light contribution to the continuum 
from the host galaxy does depend on the aperture, and it 
is necessary to correct this effect. To aco mplish this goal, 
we adopt the scheme by Peterson et al. ( 1995 ). Which is 
based on the ratio: 



Fr, 



F(4959 + 5007) 



F(4959 -I- 5007) 



G, (1) 



ohs 



where i^(4959-K5007) is the absolute flux in the [O ill] dou- 
blet and the value in brackets is the continuum to [O ill] 
lines observed flux ratio; G is an aperture dependent cor- 
rection factor to account for the host galaxy light. Since 
most of our spectra were obtained with the 6 m telescope 
through an aperture of 2"x6", this instrumental setup was 
adopted as the standard one (i.e. G = by definition). The 
value of G for the spectra obtained at the 2.1m telescope 
(aperture of 2.5" x 6.0") is G = 0.046, while the value 
obtained for the 1 m telescope (aperture of 8" x 19.8") is 
G = 0.321 ± 0.028. In Table | we fist the derived values 
of the H/3 flux in the 4990 A to 5360 A wavelength range, 
jointly with the continuum flux at 5125 AThe mean error 
(uncertainty) in our flux determination for both, the H/3 
and the continuum is « 3%. These quantities were esti- 
mated by comparing results from spectra obtained within 
a time interval shorter than 2 days. 

In Table | we list our results, there we give: 1 - The 
Julian date; 2 - a code, according to Table S 3 - F(H/3) 
the total flux (in units of 10"^'^ ergs s~^cm ); 4 - eh/S 



the H/3 flux error; 5 - Fc 
units of 10^^'^ ergs s~^cm~ 
scope aperture (2"x6"); 6 
flux error. 



the continuum flux at 5125 A(in 
^A~^), reduced to the 6 m tele- 
- Er - the estimated continuum 



3. Data analysis 

3.1. A historical B light curve 

In order to construct a fight curve in the B band dat- 
ing back to 1966, in addition to our photometric data 
for the 1997-2000 period, we collected from the literature 
all the available photometry for 3C 390.3 reported during 
the last 30 years. We included the photoelectric observa- 
tions, carried out in 1965-1967, 1971 and 1982 by Sandage 
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Table 5. Observed H/3 and continuum fluxes 
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197^; p^ ; p^8^; Selmes et al. ( |l975| ); Scott et al. 

( |1976D ; Pica et al. ( |1980D . In order to increase their accu- 
racy, the photographic observations were averaged in 5 day 
bins. We also included a few spectral observations by Perez 
et al. (1988), Lawrence et al. ( 1996| ) and the International 
Monitoring data in the framework of the AGN Watch con- 
sortium for 1994-1995 (Dietrich et al. |1998[ ) and spectral 
continuum light curve data, discussed here. The result- 
ing light curve in B band is shown in Fig. ^ There, we 
can see a noticeable increase in luminosity with signifi- 
cant fluctuations between 1965 and 1977. The amplitude 
of the outburst from the middle of 1969 until the middle 
of 1970 was about 2". In 1979 the brightness of 3C 390.3 
decreased to 17™ and remained at this level until 1983. At 



( |1973D and Neizvestny (1986); spectrophotometric obser- 
vations for 1970-1979 by Yee & Oke (|198lD, from which 



we derived B magnitudes in excellent agreement with con- 
temporary photoelectric photometry. We also make use 
of the photographic photo metry available, ie. 1967 -1980 
from Cannon et al. ( |1971 ); Babadzhanyants et al. ( |1973| ; 



that time, Pension & Perez (1984) noted that the broad 
component of H/3 had disappeared, and the spectrum of 
the object became quite similar to that of a Sy2 galaxy. 
A similar behavior was observed in NGC 4151, for which, 
after a long lasting photometric minimum between 1984 
and 1989 the broad wings of H/3 also disappeared an d was 
then classifi ed as Sy2 at that time (Pension & Perez 1984 ; 
Lyuty et al. 1984 ). In the case of 3C 390.3, after maximum 
light in 1970, on the descending branch of the light curve, 
important oscillations with an amplitude up to 1™ were 
observed. Unfortunately, we have not found photometric 
observations of this object for the 1983-1993 period and 
hence, the photometric behavior of 3C 390.3 during this 
period of time is unknown. A detailed study of the char- 
acter of the visible light variability in 3C 390.3 will be 
presented in a separate paper. 



3.2. Variability of the Hf3 emission line and the optical 
continuum 

The photometric data obtained during 1997-2000 is plot- 
ted in fig.0. The light curve shows an almost sinusoidal 
change in brightness with a maximum amplitude of about 
0™.8 in the B band. Brightness maxima ocurred in May- 
June 1998 and in March-April 2000. The light curves in 
the V and R bands are similar in shape to those of the B 
band, but have a smaller amplitude (« 0™.5), as expected 
from a blue variable continuum. In fig. |l|, one can also no- 
tice small amplitude light fiuctuations, superimposed on 
longer time-scale changes. Between 1995 and 1999 the H/3 
flux changes reached a maximum amplitude of 2.7, while 
those of the continuum at 5125 A, show a maximum am- 
plitude of « 3.2 (Table ||). These changes are evident in 
a simple inspection of the spectrum of the high-activity 
state (20 Mar 1996) as compared with corresponding one 
of the low-activity state (9 Sep 1997). These are plotted 
in Figure ^. From the combined light curve, (Figj^). For 
the 1995-1999 time interval, we confirm a maximum am- 
plitude flux variations in B of « 3.2. This value is in excel- 
lent agreement with the value derived from spectroscopic 
continuum fluxes, rendering further support to our basic 
assumption: ie. that the flux of the [O ill] lines did not 



8 



Shapovalova, A.I., Burenkov, A.N., Carrasco, L., et al.: H/3 monitoring of 3C 390.3 



3^ 

« 2^ 
u. 

3^ 
>< 2^ 



1 

5 



- a) 



- b) 



1^1 - 



H 5 



J 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 I 1 L 



1 ^aS- 



I.I.I .1.1 . I . I . I .1.1 



t 



30000 30300 31000 31300 

JD -2.100000 



I 



6 

o 

m 
w 

Sh 

0) 



CO 

T— I 

I 



1^ 



1 1 1 1 1 1 1 1 1 1 1 1 1 

- (a) 0°^ - 

- °d> - 

o o 

- - 

; °o r = 0.67 : 

~i 1 1 1 1 1 1 1 1 r 


1 1 1 1 1 1 1 1 1 1 1 1 1 

-(b) 8 - 

- o° o y° - 

o o 

- ir °° - 

: f °e : 

: o r = 0.74 : 

~i 1 1 1 1 1 1 1 1 r 


-(c) °o. 

- o ° - 
o y 

o / o 

o ° 

r = 0.82 - 

1 1 1 1 1 1 1 1 1 1 1 1 1 


: ('i) - 

o ^ o 

o r = 0.94 " 

1 1 1 1 1 1 1 1 1 1 1 1 1 



, , -15 -2 -1 -1, 

F(5125 A) (10 ergs cm s A ) 

Fig. 5. The F(H^) and Fc(5125) plotted against each 
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gressions are also given. 



Fig. 4. 3C 390.3 Light curves of the H/3 panel a and the 
5125A optical continuum emission panel b during 1995- 
2000. Filled circles represent the spectral data. The pho- 
tometric points are represented by open circles. Units 
are 10~^^ ergs cm~^s~^ for H/3 fluxes and 10"'" 



ergs 



cm^^s^^A^^) for the continuum ones. In panel c we plot 
- Vred - Vbluc, the vclocity difference of the peak emission 
of the red and blue bumps (in km s~^), derived from the 
H/3 difference profiles 



change during the period of time covered by our observa- 
tions. 

Figures ^ and ^3 show the light curves for the inte- 
grated H/3 (4990-5360 A ) and continuum (5125 A) fluxes, 
these data are given in Table 5. In order to improve the 
time resolution of our continuum data set, we have added 
some photometric points in the V band to our spectral 
continuum data. The V band data was converted into 
flux adopting the calibration by Johnson ( 1966D . A com- 
parison of the spectral continuum photometry at 5125 A 
(Fc(5125)) with 10 simultaneous observations in the V 
band, yields the following transformation equation: 



i^c(5125) ^ (0.941 ± 0.139) x Fx{Y) 
(correlation coefficient r 



(0.722 ±0.326). (2) 



given in units of 10 ergs cm ^s 



0.923). Here the fluxes are 
~'A-')- 

In Figures Qa and ||b two outbursts in H/3 and the con- 
tinuum are evident, with time intervals between minima 
of « 1000 days (Oct. 1994 - Jul. 1997) and w 700 days 
(Jul. 1997 - Jun. 1999), respectively. Also, the presence of 



a time lag between the continuum light changes and the 
response of the H/3 line, is also noticeable. 

Fig. ^,b,c,d shows the observed correlation between 
the continuum and H/3 light curves, for a time lag between 
the continuum and line variations of 0, 20, 40 and 100 
days respectively. There, it is clearly seen that the variance 
of the data is minimal while the correlation coefficient is 
maximal (r=0.94) for a delay of 100 days. In latter case, 
wc obtain the following relationship: 

F(H/3) = (0.957 ± 0.056) x i^c(5125) -I- (0.640 ± 0.102) (3) 



Here i^c(5125) fluxes are expressed in units of 
10~^^ ergs cm~^ s~^ A~^ and the F(H/3) fluxes in units 
of 10"^'^ ergs cm^^ s^^. The values chosen for the delays 
are typical values previuosly mentioned in the literature. 
In Sect. 3.4 we present a detailed discussion of the delays 
from a cross-correlation analysis of our data. 

3.3. The H/3 wings 

In order to investigate the flux in the different parts of the 
H/3 profile, we studied three different velocity bins: a blue 
wing (-7100km s^^ to -1900km s^^), a core (-1900 
km to -1-1900 km s^^) and a red wing (-1-1900 km s~^ 
to -1-7100 km s~^). The light curves for these are presented 
in Figures ^ a,b,c, there, we can see that the changes in 
both the blue and red wings and the core occur quasi- 
siniultaneously. This fact implies that, there are predom- 
inantly circular motions in the region where the broad 
H/3 emission originates. The blue-to-red wing flux ratio 
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Fig. 6. The light curves, (a)-(c) for the blue wing in the 
5010-5099 A (-7100 km s^^ to -1900 km s^i), interval, 
for the core in the 5099-5164 A (-1900 km s'^ to +1900 
km s"'^), interval, and for the red wing of H/3 emission line 
in the 5164-5253 A (+1900 km s^^ to +7100 km s^^), in- 
terval, respectively. Panel (d) shows the temporal behavior 
of the the blue-to-red wing flux ratio R 



R is also shown in Fig. |^. Veilleux & Zheng ( |l99l| ) no- 
ticed that between 1974 and 1988, the blue-to-red wing 
flux ratio followed a fairly smooth, almost sinusoidal pat- 
tern with an apparent period of 10.4 yrs. This trend was 
confirmed by Bochkarev et al. (1997a) with observations 



that extended until 1995. Although maxima for the ratio 
R were observed in the years: 1975, 1985, and 1995, the 
sinusoidal trend of R did not continue during the 1996- 
1999 period. In Fig. ||d, it is evident that R continued to 
increase in a nearly monotonic fashion between 1995 and 
1999. 



3.4. Cross-Correlation Analysis 

3.4.1. The CCF of the H/3 broad and continuum 
emissions 

In order to determine a more accurate time lag than the 
one described in 3.2, we carried out a time series analysis 
using the cross-correlation function (CCF). The CCF was 



calculated by means of the interpolation method described 
by Gaskell & Sparke ( |1986D and White & Peterson (cite- 
whi) . We computed both the lags related to the CCF peak 
(tmax) and the CCF centroid (teen)- The value of {teen) is 
that of the center of mass of the CCF for positive delay 
values. The CCF is calculated by pairing each of the real 
data points in both time-series with linearly interpolated 
points for any arbitrary time delay. 



According to Gaskell and Peterson (1987) not every time 
series is suitable for the cross-correlation analysis, only 
those series, for which the autocorrelation function (ACF) 
width is wider (by at least 10%) than the width of the cor- 
responding function for the sampling window (ACFW), 
contain relevant information about time lags. The ACFW 
is computed by repeated sampling of white noise light 
curves, convolved with the observational window. This 
procedure provides a measure of how much of the width of 
the ACF is due to the interpolating scheme, instead of the 
real correlation width of a continuous time series sampled 
at discrete times. In addition, the half width of the ACF at 
a zero correlation level determines the characteristic size 
of the region from which variations originate, setting an 
upper limit to the BLR size. 

The ACFW was calculated for the time series of the obser- 
vations, upon which, random Gaussian noise correspond- 
ing to the observing errors in the time series was superim- 
posed. The mean ACFW for one hundred of those realiza- 
tions is shown in Fig. There we can see that the ACFW 
is much narrower than the ACF for either the continuum 
or the up emissions. 

Usually, the time delay is determined either by the lo- 
cation of the peak, or that of the center of mass of the 
cross-correlation function between the continuum and the 



emission line light curves. Robinson and Perez ( 1990 ) ar- 
gued that the position of the CCF peak {tmax) yields the 
time delay corresponding to the inner radius of the BLR, 
while the position of the center of mass of the CCF {teen) 
relates to the size of the BLR weighted by luminosity. 
The errors, in the determination of the time lags, were 
estimated through a Monte Carlo simulation. From 1000 
independent realizations of the CCF, the tmax and teen 
distributions were obtained. Prior to the computation, the 
time-series involved, were corrupted with Gaussian noise 
in correspondence to the typical errors in them. Cutting 
the distribution functions at a 17% probability level, the 
uncertainty related to a 67% confidence level {la error) 
was determined. 

Given that, during the International Monitoring 
Programme of the Sy galaxy NGC 5548, for a pe- 
riod of 8 years, it was observed by Peterson et al. (1999) 



that the time lag in this object does not have a constant 
value, and that changes of the delay may be related 
to the luminosity variations of the central source. And 
since, during time interval cover by our observations, 
the photometric and spectral light curves of 3C 390.3 
show two distinct fiares. It is important to determine the 
delays not only for the entire time-series, but also for the 
segments of it, that correspond to relevant events in the 
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Fig. 7. The Auto- and Cross-Correlation Functions. In panel (a), the narrow eurve represents the sampling window 
autocorrelation function (ACFW) for the spectral continuum and H/3 line. Also, the autocorrelation functions (ACF) 
for the H/3 label (1) - and the continuum light curves (2) are shown. On panels (b) and (d), the CCFs for the periods 
1995-1999 and 1997-1999 respectively are shown. Curves label (1) represent the CCF between the H/3 flux and the 
spectral continuum at 5125A; curves labeled (2) represent the CCF for the H/3 flux and the combined continuum and 
label (3) refers to the CCF computed as in case (2), but for a restricted interpolation of 100 days. Panel (c) shows 
the CCF between: (1) - H/3 blue and red wings light curves , (2) - H/3 red wing and H/3 core fluxes and (3) - H/3 
blue wing and H/3 core light curves. In every panel the units are: delay in days (abscissae) and correlation coefficient 
(ordinates). 



light curves. Therefore, analysis of the light curves was 
carried out for various cases, as follows: 

The CCF analysis has been carried out separately for 
the entire data set - case A (JD=2449832-2451527), and 
for the two sub-sets (case B, JD=2449832-2450701; and 
case C JD=2450618-2451527). It should be mentioned 



that for case B, we do not have supplementary photom- 
etry. While in case C we have included our broad band 
photometry. For this reason, the flux cross-correlation was 
calculated not only for the spectral contimmm fluxes Fgc, 
but also for the combined continuum data Fscv To the 
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latter, we added our photometric data in the V band to 
the spectral continuum data as discussed above in 3.2. 

Since many spectral observations were obtained with 
long time intervals between them, the CCF was calculated 
restricting the interpolation processes to time intervals no 
longer than 100 days (cases - AD, CD). That is, prior to 
the CCF computation, both time series were analyzed as 
to whether or not they were suitable for pairing. With this 
100 day restriction in the pairing, possible effects associ- 
ated to long time gaps are reduced. 

The results of the cross-correlation analysis are pre- 
sented in fig. 1^ and summarized in Table ^. There we 
list: in Col. 1 the time-series case; in col. 2 - the cross- 
correlation case: (F(H/3) - Fgc) - between the continuum 
flux F(5125) and H/3 flux; (F(H/3) - Fscv) - between the 
combined continuum flux (F(5125)-|-F(V)) and H/3 flux; 
in col. 3 ~ tmax ~ the time lag in days, determined from 
the position of the CCF maxima; in col. A - la interval for 
tmax] in col. 5 - rmax~ the cross-correlation coefficient for 
the CCF maxima; in col. 6 -teen ~ the time lag in days, 
determined from the centroids of the CCF; in col. 7 - the 
estimated la uncertainty for teen] in cols. 8 and 9 - the 
number of points, used in the computations - N(H/3) and 
-Nc. 

Inspection of Fig. ^ and Table ^, shows that: 

1. All cases yield a time lag of H/3 relative the continuum 
flux changes tmax ~ 100 ± 10 days or a little less for 
teen « 90 ± 10 days. 

2. In the time interval from Jun.1997 - Dec. 1999 (case C), 
an additional, shorter delay of 33-37 days is apparent. 

Our delay values differ from the results of Dietrich et 
al.(1998). From Fig. | one can clearly see, that the 20 
day lag, obtained by these authors yields a much larger 
variance and lower correlation coefficient between line and 
continuum fluxes than the 100 day lag. In our data we do 
not see a lag of 20 days or something like it. The reason 
of this discrepancy is not clear to us. Perhaps it is related 
with the longer time interval covered by our observations. 
It is conceivable that the time delay in this object does not 
have a constant value, and that changes of the time lag 
may be related to the luminosity variations of the central 
source. A similar case to th at of the Sy galaxy NGC 5548 
studied by Peterson et al. ( 1999| ). 



3.4.2. The CCF for different parts of the H/3 profile 

We calculated the cross-correlation of the responses of the 
blue wing (-7100 to -1900 km/s), of the core (-1900 to 
-hl900 km/s) and of the red wing (+1900 to +7100 km/s) 
of the H/3 profile to variations of the optical continuum. 
The delays found were 100-113 light days for tmax, and 
86-102 light days for teen- However, it should be men- 
tioned that 7% of the maxima in the distribution of the 
CCF peak for the case of the F(H/3) blue wing - Fgev cor- 
respond to a delay of about 40 days, almost the same as 
in case C previously discussed (see Sect. 3.4.1). 




-10000-5000 5000 10000 
Velocity [km/s] 

Fig. 9. The mean and rms H/3 profiles in velocity units 



The possible differences in the response of the various 
parts of the line profile can be studied through their cross- 
correlation functions. In our case, the cross-correlation 
of the blue-core, red-core and blue-red H/3 components 
are presented in Fig. The corresponding time lags for 
the blue and red wings relative to the core, and of the 
wings relative to each other are listed in Table |[ There, 
we can see that our CCF analysis for the H/3 wings does 
not revealed any delay in the variations of the line wings 
with respect to the central part of the line, or relative to 
each other. These results preclude the possibility of having 
a BLR velocity field dominated by radial motions. 



3.5. Hf3 profile changes 

In Fig. H we present the profiles of the broad H/3 compo- 
nent, in velocity units relative to the narrow component. 
There the narrow [O ill] and H/3 emissions were removed 
by means of the Gaussian-fit procedure previously de- 
scribed. Characteristic features in the profile of the broad 
H/3 line are the presence of a blue bump and an extended 
red wing. The blue wing of H/3 was brighter than the red 
one during the 1995-1999 interval. In 1996, the intensity 
of the wings was very strong. Yet in 1997, the emission in 
the H/3 wings was rather weak. 
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Fig. 8. Typical broad H/3 profiles (solid line) after continuum and narrow components subtraction (dashed line). 

3.5.1. Mean and Root-Mean-Square Spectra 1995-1999 lapse. From Fig. | (bottom) the FWHM-(rms) 

is about 12 000 km s"^ 



The comparison between the average and root-mean- 
square (rms) spectra provides us with a good measure of 
profile variability. Since, any constant contribution to the 
spectra is excluded from the rms spectrum. The mean H/3 
profile after removing narrow the H/3 and [O ill] lines and 
the absolute rms variations per unit wavelength are shown 
in Fig. ^ Our results are similar to those of Dietrich et 
al. (1998): ie. the mean H/3 profile is clearly asymmet- 
rical with a full width at zero intensity of about 20 000 
km s~^. The blue bump in the mean spectrum is located 
between -3000 and -4000 km s^^ while the red bump is 
seen between -1-4000 and -1-5000 km s^^. The mean profile 
shows a blue bump brighter than the red one during the 



3.6. The velocities of the Hp blue bump 

We have measured the radial velocities of the emission 
peak of the H/3 blue bump relative to the narrow H/3 com- 
ponent from Gaussian fits of the top of the blue bump 
down to the 30% level of the peak brightness. A method 
similar to that applied by Eracleous et al. (1997). A com- 
parison of our blue bump radial velocities with theirs, for 
the same epochs show a very good agreement between 
both sets of observations. The mean difference being about 
±60km s~^, a value close to the uncertainties of the proce- 
dures adopted. Our results are presented, in graphic form 
(open circles) , in Figure |l^. Each point represents an av- 
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Time 


Case 




la 




teen 


la 


N(H/3) 


Nc 


series 






interv. 






interv. 






A 


F(H/?)-Fsc 


105 




0.928 


99 




44 


44 


A 


F(H/3)-Fscv 


100 


97—101 


0.900 


86 


79—92 


44 


127 


AD 


F(H/3)-Fecv 


99 




0.930 


98 




37 


45 


B 


F(H/3)-Fsc 


98 




0.946 


109 




20 


21 


B 


F(H/3)-Fscv 


99 


93—113 


0.954 


112 


98—114 


20 


29 


C 


F(H/3)-Fec 


105 




0.885 


90 




27 


26 






34 




0.800 










C 


F(H/3)-Fecv 


34 


33—37 


0.829 


63 


50—72 


27 


109 






99 


97—119 


0.784 










CD 


F(H/3)-Fecv 


99 




0.859 


78 




24 


33 






33 




0.803 










A 


F(H/3)b-F(H/3)c 





-3.7—1.2 


0.809 


-1.6 


-5.0—1.2 


44 


44 


A 


F(H/3)r-F(H/?)c 


-1.0 


-3.3—0.9 


0.957 


-1.5 


-4.6—0.8 


44 


44 


A 


F(H/3)b-F(H/3)r 


0.5 


-1.3—7.8 


0.894 


0.6 


-2.1—8.0 


44 


44 



erage of several spectra obtained during a lapse of a few 
months. The absolute value of the radial velocity of the 
blue bump presents a minimum, observed in 1997, at « 
-2900 kms-\ and a maximum occuring in October 1999 
at ~ —3700 km s^^ . The velocity shift of the blue bump 
reached a value of about 800 km during this time. 
The changes of the blue bump velocity between 1967 and 
1999, taken from the literature, are also presented in Fig. 
[lO| . The short dashed line represents the best fit of the 
double- line spectroscopic binary model by Gaskell ( 1996D , 
while, the long dashed line represents the fit by Eracleous 



et al. (1997). Onecan see that our data are reasonably well 
described by the latter. These authors have shown that, 
if the displacement of the blue bump peak velocity origi- 
nated from individual broad-line regions, associated with a 
massive binary BH, then the infered rotation period would 
be about 800 years, and the associated mass would then 
be larger than IO^^Mq. This value for the BH mass is 
much larger than the expected by comparison with other 
BH candidates found at the center of galaxies. Therefore, 
these authors rejected the binary BH interpretation. Our 
results further support their conclusions. 



-3000 



-3500 - 



-4000 - 



-4500 - 



-5000 - 




1965 



Fig. 10. The average radial velocity curve of the blue 
bump peak. The different symbols represent: triangles - 
the annual H/3 average data from Gaskell (199£), open 
squares - the Ha data from Eracleous et al. (1997), open 
circles - our H/3 data. The different best fit models are 
explained in the top left panel. 



3. 7. The H/3 difference profile 

We obtained th H/3 difference profiles by subtraction of the 
minimum-activity average spectrum (9 September, 1997) 
from individual spectra. This is done after continuum, nar- 
row H/3 and [O ill] emission had been substracted as well. 
When the spectra for subsequent nights, presented small 
differences ie. H/3 integral flux remained constant within 
3 to 5%, they were averaged and the spectra of low signal 
to noise ratio (S/N< 20) were excluded from the analy- 
sis. The annual average differential profiles of the broad 
H/3 line are shown in fig. |l^. In order to determine the 
blue and red peak velocities, the mean difference profiles 
were fitted with three Gaussian functions; one to the blue 
bump, one to the core, and one to the red bump. The de- 
rived differences, of the peak velocity for the red and blue 



bumps, Vr-Vb, are plotted in Fig. ^ c. Despite the errors 
in the determination of the velocities, a distinct anticorre- 
lation is observed between the difference Vr - Vb and the 
fiux changes in both H/3 and continuum emission, as it is 
clearly seen in Fig. ^. 

The annual mean velocity of the blue, red and core 
components derived from the Gaussian analysis, jointly 
with the Vr-Vb velocity difference of the bumps, are listed 
in Tabic 0. There we notice that, the radial velocity of the 
blue component increased in absolute value from —3200 
km s^i in 1995-1996 to -5200 km s^^ in 1999. At the 
same time the radial velocity of the red component in- 
creased from -f-4900 km s'^ in 1995-1996 to -1-7000 km s^^ 
in 1999. Then, the absolute velocity of both components 
increased by about 2000 km s^^ during this period of 
time. That is, the difference Vi-Vb increased by about 



14 



Shapovalova, A.I., Burenkov, A.N., Carrasco, L., et al.: H/3 monitoring of 3C 390.3 



Table 7. Radial velocities of the blue bump, red bump and core of H/? derived from the difference profiles, jointly 
with contemporary H/3 fluxes. 



Velocity (km s ^) 




Time Interval 






Fluxes 


Apr. 1995-Nov. 1996 


(Mar .-Aug.) 1997 


(Jan.-Oct)1998 


(Jul.-Oct.)1999 


Vbluc 


-3204±187 




-3867±75 


-5215±118 


Vrcd 


4899±155 




5361±134 


7001±106 


Vcorc 


617±187 


76±249 


335±235 


1170±244 


Vr-Vb 


8103±243 




9228±154 


12216±158 


F{m* 


2.804±0.37 


1.733±0.14 


2.367±0.21 


1.932±0.25 


con 


2.132±0.54 


1.035±0.13 


1.611±0.14 


1.406±0.17 



* - in units of 10 ^'^ ergs s ^cm ^ 
** - in units of 10~^^ ergs s~^cm~^A 



4000 km s~^. Table ^ lists the observed mean fluxes of 
H/3 (F(H/3)) and continuum {Fcon) light in the same time 
intervals. It is evident that the lower annual average ve- 
locities of the blue and red components, or their difference 
correspond to larger mean flux values. In the frame of an 
accretion disk model, the optical continuum changes are 
a consequence of changes of the luminosity of the central 
X-ray source, causing a variable disk irradiation. 

The observed velocity variations could be explained if, 
at different times, the zone that contributes with the max- 
imum luminosity to H/3, changes in radius within the disk. 
The position of this zone will depend on the luminosity of 
the central source. Then, it is likely that during 1995-96, 
we were observing enhanced H/3 emission from parts of the 
disk which are located further away from the BH ie, when 
the luminosity of the central source is higher, and conse- 
quently heated the outer regions of the disk. This results 
in a maximum emission zone located at larger disk radii. 
Later on, in 1999, maximum emission came from parts of 
the disk located at a smaller radii, when the luminosity of 
the central source was much lower. In the case of a circular 
disk the ratio of the maximum H/3 emission radii must be 
proportional to the inverse squared velocity. Hence, if we 
consider the velocity values listed in Table |^ for 1995-96 
and the average velocity values for 1998-99, the ratio of 
the radii of maximum H/3 emission is about 1.77 ± 0.2. 
This is in excellent agreement with the ratio of the size of 
the emitting zones as inferred from the lag time ratio of 
the centroids of the CCF (ie. 1.78 ± 0.3, see Table |). 

These transient phenomena are expected to result from 
the variable accretion rate close to the central black hole. 
Note that from the theory of thin accretion disks (Shakura 
& Sunyaev 1973) one can calculate the energy dissipated 
through viscosity, the resulting spectrum emitted by the 
disk, and time scales in which, the emission from the disk 
varies in response to a variable accretion rate. These are 
rather long (sa 10^ years for a 10 ""Mq BH ). However, 
as the accreted matter approaches the BH, at a distance 
of few tens of gravitational radii, a large fraction of the 
accretion energy is transformed into X-ray radiation. The 
rapid variability of this radiation on time scales of hours. 



days and weeks is due to changes in the final accretion 
rates. 

Currently, a model in which the surface temperature of 
the disk is modulated by the irradiation caused by a vari- 
able luminosity of the X-ray source has been suggested by 



several authors (see Ulrich, 2000 and references therein). 
In some models it is assumed that the rapid variability 
(days, weeks) of the continuum X-ray source can be ex- 
plained by the fact that a significant fraction of the ac- 
cretion energy near the black hole is spent on heating a 
corona by a mechanism reminiscent of flares in the solar 
corona (Galeev et al. 1979). Hot electrons, in turn, transfer 
some of their energy to ambient soft X-ray and UV pho- 
tons (some emitted by the disk), producing medium en- 
ergy and hard X-ray radiation via inverse Compton emis- 
sion. 



3.8. Modeling of the H(3 Broad Components 

As mentioned above, several models have been proposed to 
explain the double-peaked line profiles. Among them, the 
emission by an accretion disk or a torus, whose presence 
seems essential to fuel a black hole. Our study favours the 
formation of the broad H/3 line of 3C 390.3 in an accretion 
disk. Superluminal motions have been observed at radio 



frequencies (Alef et al. 1996), these are indicative of the 
presence of a relativistic jet, with its axis aligned close 
to the line of sight (Orr & Browne, 1982). Therefore, the 
disk's axis is probably close to the line of sight. 

Under the assumption of an accretion disk as the site 
of emission line formation, the study of the profile changes 
with time, may provide relevant information, even if we do 
not advance any hypothesis about the disk structure or to 
how it is illuminated. Moreover, we could reject the disk 
hypothesis, if we find some inconsistency in the results (for 
instance a change of the BH mass or the disk inclination) . 

We assumed that: i) the Balmer lines are emitted by 
a disk rotating about a central object with Keplerian ve- 
locity, and ii) the disk emissivity varies with the distance 
r = R/Rg to the center, where Rg is the gravitational ra- 
dius {2GM/&') with M being the mass of the BH. Then we 
computed the profiles emitted by such disk, taking into ac- 
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Fig. 11. Annual averages from 1995 to 1999 of the differ- 
ence profiles for the broad H/3 line. 



count the relativistic Doppler and gravitational redshifts, 
using the eq uation s given by Gerbal & Pelat ( 1981 ) and 
Chen et al. ( 1989 ). These equations are valid for a disk 
seen at inclination angles i smaller than 80deg, which is 
the case for 3C390.3. The resulting profiles are double- 
peaked and asymmetric. They were convolved with the 
instrumental profiles and fitted to the observed ones, us- 
ing the Davidon- Fletcher-Powell method found in Minuit's 
package (James 1994). 

One important parameter in the fitting procedure is 
the inclination angle i, which must remain constant. For a 
simple power-law for the disk emissivity ie. F{Il(3) oc r"'', 
the line wings cannot be properly reproduced, a fact al- 
ready noticed by Chen et al. (1989) for Arp 102. The pro- 
file fits are improved significantly, when a double power- 
law emissivity function is adopted. In this case, the im- 
portant fitting parameters of the model are: the indices qi 
and q2, the radius ri at which the slope changes from qi 
to (?2j and the outer radius of the disk, rout (all radii being 
expressed in units oi Rg). The inner disk radius does not 
play an important role. The four relevant parameters are 
expected to vary with the continuum flux level and the 
illumination of the disk. 

In this paper, we present only a few of the results, a 
more detailed description of the model and the results for 
the Ha and H/3 profiles will be presented in a separate 
paper. 



The fits of four selected H/3 profiles, distributed in 
time, and representative of maximum and minimum ac- 
tivity states, are shown in Fig. ^ From our fitting pro- 
cedure we find that: i) z = 25 ± 3 degrees, consistent 
with a constant value of i, and ii) surprisingly, despite 
the large changes in the continuum and line fluxes, qi ~ 1 
and ri « 200i?g do not vary significantly, while 52 varies 
between 3.46 and 3.51 for the fits to the four profiles. 

One should seek a physical explanation of this behav- 
ior. According to the results, the bulk of the line emission 
is produced at a radius close to ri i.e. at i? « 200i?g. 

It is interesting to note that this value corresponds to 
the radius where the accretion disk b ecomes gravitation- 
ally unstable (cf. Collin & Hure 2001 ). At larger radii, the 
emissivity decreases very rapidly, faster than it would be 
expected in the case a continuum source illuminating the 
accretion disk, or a system of rotating clouds having a 
constant density with a canonical BLR value of n > 10^ 
cm^'^. In this case, i^(H/3) would roughly be oc (cf. 



for instance Dumont, et al. 1998), where Fine is the flux 



incident on the disk. And since at large distances from 



the continuum source Fine oc 



one would find that 



-1.5 



for a constant density. If the disk is geo- 



metrically thick (a "torus"), the exponent is smaller. Fine 
being proportional roughly to r~^. This means that the 
density decreases rapidly with increasing radius, as about 
r~^, as expected if the disk is in a state of marginal in- 
stability (Collin & Hure, 1999). On the other hand, for 
smaller radii, the emissivity varies slowly with r, which is 
compatible with the disk irradiated by the central source 
of continuum if the density is larger than 10^ cm^'^. Thus 
this simple change of the density in the accretion disk 
could induce a "physical radius" of the emission region 
quite independent of the value of the continuum flux, and 
could explain why the radius of maximum emission does 
not vary strongly during the monitoring. 

So, our study favours the formation of the broad H/3 
line of 3C 390.3 in an accretion disk. 

3.9. The mass of the Black Hole in 3C 390.3 

Using the virial theorem jointly with the reverberation- 



rms method (cf. Peterson et al. 1998), the mass of the 
black hole in 3C390.3 can be estimated from the relation- 
ship of Wandel et al., (|l99|): 



Mr 



(1.45 X IO^Mq) 



CT 



It — day 



rms,3 ' 



(4) 



where Wrms,3 = w_FWHM(rms)/10^ kms"^ 

Using our results: ypwHAii'^'T^s) ~ 12 000kms~-^ 
(Fig. |b), (cr/lt-day) « 100 (Table 0), we obtained 
Mj-ev ~ 2.1 X 10^ Mq. This value is substantially larger 
than M„p„ 



s 4 X 10** Mq, obtained by Wandel et al. 
( |1999|) . The difference is due to the fact that the deter- 
mination of Wandel et al. is based upon: 1) a time lag of 
24 days, instead of the 100-day lag found in the present 
study; 2) a the value for vfhwhm (njis) ~ 10500 km s~^ 
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a) ^21.0j,1996) 




Fig. 12. The fits to the H/3 observed profiles for differ- 
ent epochs, by an accretion disk model with an inclina- 
tion angle of about 25°, and with a region of maximum 
emission at about 200 gravitation radius. The dotted line 
represents the model profile, while the full thin line the ob- 
served one. In each window, the lower panel represents the 
residuals for the following dates: JD=24410163, 24410701, 
24411025, 24411426. 



instead of our value for VFWHMirias) « 12 000kms ^. 
The upper limit for the mass of the BH in 3C 390.3, esti- 
mated from X-ra y var iability is M{BH) < 1.5 x 10^° Mq 
(Eracleous et al. |1997[ ). 

s M{BH), we estimate a value 



Assuming that Mre 
for Rg (the BH gravitational radius), for the derived value 
of the BH mass: 

Rg = 2GM{BH)/c^ « 6.2 x lO^-^cm. 

According to our model calculations (Sect. 3.8 ), the 
bulk of H/3 emission is produced at r w 200 Rg, which 
corresponds to a time lag r = r/c « 48 light days. This 
implies a shorter time lag than the one derived in this 
paper (Table Can these values be reconciled? The re- 
verberation method assumes that the FWHM of the lines 
reflects virialized motions, while we assumed that the mo- 
tions are purely rotational. These two working hypotheses 
may explain differences in the size of the emitting region 
on the order of w 40 - 50 %. Then the time lag for 200 Rg 
could be about r « 70 It. days. A value which is closer to 
the mean lags listed in Table 0: Teen ~ (86 ± 24); and 
Tmax ~ (67 ± 32). Worth, mentioning is the good agree- 



ment (within 20 "cool dense material" reprocesses X-ray 
radiation into the Fe Kalpha line detected by Eracleous, 
et al. (|199|). 



4. Conclusions 

The results of a 4-year (1995-1999) spectroscopic and 
broad-band BVRI photometric monitoring of the AGN 
3C 390.3 are presented in this paper for the H/3 spectral 
region. 

Our main conclusions are: 

1. The historical light curve in B shows a large increase of 
brightness during 1970-1971 followed by a gradual de- 
crease until 1982. At that time, the nucleus of 3C 390.3 
was in a minimum luminosity state. On the descending 
branch of the light curve, flares with an amplitude up 
to 1™ were observed. 

2. The broad H/3 and continuum fluxes vary by a factor 
of about three during the time interval 1995-1999. 

3. Two large amplitude outbursts in the intensity of the 
H/3 line and the continuum, with different duration 
were observed: (Fig. |): in (Oct. 1994- Jul.1997) - a 
brighter flare lasting wlOOO days and another one in 
(Jul. 1997 - Jun. 1999) with that lasted «700 days. 
The event duration is defined as the interval between 
minima in the light curve. 

4. From a cross-correlation analysis, a delay of « 100 days 
in the response of the H/3 emission with respect to the 
optical continuum changes is obtained from the obser- 
vations at any given time. Yet, for observations ob- 
tained between the June 1997 until 1999 (JD2450618 - 
2451527 - case C), a period in which the sampling was 
better, a shorter lag of 33 to 36 days appears jointly 
with the one of « 100 days. In the light of these re- 
sults, we have to admit that either, our data set is not 
suitable for an accurate determination of the delays 
due to a somewhat poor temporal sampling, or that 
the delay for 3C 390.3 varies with time, i.e. the cor- 
relation between the flux variations in the continuum 
and in H/3 at different times is not the same. In the 
8 year mon itoring campaign of NGC 5548, Peterson 
et al. (1999) have determined, that the time delay had 
varied. They suggested that changes in the time delays 
may be related to luminosity variations of the central 
source. Yet, in other investigation of NGC 5548 (Netzer 
and Maoz, 1990) have shown, that continuum varia- 
tions with a long time scale -long lasting flares- yield 
longer time lags. Both posibilities are present in our 
results for 3C 390.3 ie.: a) the continuum outburst of 
1994-1997 with larger amplitude and duration (w 1000 
days) yields a longer time lag (100 days); b) the con- 
tinuum outburst in 1997-1999 of a smaller amplitude 
and duration (« 700 days) yields two lags: 100 and 
35 days. We believe that both, the optical continuum 
and H/3 flux variations reflect the changes in the X-ray 
irradiation, modulated by accretion rate changes. And 
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that the radius H/3 emission zone, changes due surface 
temperature changes in the disk (Uhich, ^000 ). 
The observed flux for the blue and red wings of the 
profile varied quasi-simultaneously (Fig. ^ (i.e. when- 
ever the flux in the blue wing increases, that in the 
red wing also increases). We have not detected a de- 
lay in the changes of the H/3 blue wing (—7100km s^^ 
to -1900km s-i) and of the red wing (-f 1900km s^^ 
to -1-7100 km s~^) with respect to the core of the line 
(-1900km s~^ to -1-1900 km s"^) or relative to each 
othe r. A similar result was obtained by Dietrich et 
al. ( 1998 ) from monitoring data for 3C 390.3 during 
1994-1995. These results exclude the presence of con- 
siderable radial motions in the BLR in this object. 
The blue-to-red wing flux ratio shows a monotonic in- 
crease with time for the period covered by our observa- 
tions (Fig. |d). The quasi-sinusoidal changes (period 
« 10 years) of the b lue-to -red bump flux ratio detected 
by Veilleux et al. (1991) are not present in our data 
for the 1995-1999 lapse. 

The shape of the H/3 profile has changed drastically: 
during 1996, the blue and red bump emission was very 
strong, while when the source was faintest (1997) the 
profiles had rather fiat tops (Fig. ||) . The blue wing was 
brighter than the red one during the observed period. 
The changes of the radial velocity of the blue bump rel- 
ative to the H/3 narrow emission in 1995 until Oct. 1999 



follow the same trend found by Eracleous et al. (1997). 
Their model for a binary BH with a very large mass 
value (> IO^^Mq) seems to fit our data as well. These 
authors argued that such a large value for the binary 
BH mass, is difficult to reconcile with other observa- 
tions and theory. Therefore, they reject the hypothesis 
of a binary black hole for 3C 390.3. Our results pro- 
vide further support to the dismissal of the binary BH 
hypothesis for 3C 390.3, on the basis of the masses 
required. 

The difference profiles of the H/3 broad emission re- 
vealed that the radial velocities of the blue and red 
bumps in the profiles varied, showing an increase in 
absolute value by about 2000 km s^^ during the ob- 
served period. From these results (Table 0, Fig. pT| ), 
it is infered that an annular zone in the accretion 
disk with enhanced H/3 emission could have moved to- 
wards smaller radii between 1995 and 1999, i.e. larger 
Keplerian velocities as the fiux in the continuum de- 
creased (Table ^ . These results are in agreement with 
a model in which the surface temperature of the disk 
is modulated by the irradiation of a variable X-ray 
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source (see Ulrich, 2000 and references threin). The 
X-ray variability on time scales of hours, days, weeks 
is, in turn, modulated by accretion rate changes. Thus, 
the peak velocity shifts as different profiles are also ex- 
pected to result from the variable accretion rate near 
the black hole. 

We have fitted the broad H/3 line profiles in the 
framework of an accretion disk model. We obtained 
a satisfactory agreement for a disk with an inclina- 
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tion angle of about 25° and with the region of max- 
imum emission located at about 200 Rg. The radius 
of this zone is comparable to the one in which the Fe 
Kalpha line emission originates according to Eracleous 
et al. ( |1996D . 

The mass of the black hole in 3C390.3 estimated by 
the combined reverberation-rms method (Wandel, et 
al., 119991) is AUev « 2.1 x 10^ Mq. 



Our results do not support either the models of out- 
flowing biconical gas streams or those of supermassive bi- 
nary BHs. Instead, we conclude that our results favor the 
accretion disk model. 

The velocities of the blue and red bumps and their dif- 
ference, obtained from H/3 difference proflles, show a dis- 
tinct anticorrelation with continuum flux changes. Taking 
into consideration the time lag and the H/3 fiux, we find 
that the zone of maximal contribution to line emission, 
moves across the face of the disk. It is located at smaller 
radii when the flux in the continuum decreases (bump 
velocities increase), and to larger radii when the contin- 
uum flux increases (bump velocities decrease) . These tran- 
sient phenomena are expected to result from the variable 
rate of accretion close to the black hole. At low accretion 
rates, the flux from the central source may decrease by 
a large factor, and the integral flux in the line could de- 
crease. Then, the maximum emission zone of the broad 
Balmer lines could shrink to such small radii, that the 
lines could become extremely broad and of low-contrast. 
Consequently, at those times, they would be undetactable. 
This may be the reason why the broad emission lines were 
not s een in the 3C 390.3 spectrum in 1980 (Heckma n et a l. 
1981 ) or were very weak in 1984 (Penston & Perez 1984 ). 
At those times the central source brightness was at long 
lasting minima (Fig. |^). 

It is intresting to note, that the model calculations of 
Nicastro (120001 ) point to a relationship between FWHM 
broad emission lines and accretion rate for a given BH 
mass — the lower the accretion rate the greater the line 
width. For different BH masses, there is a minimum value 
of the accretion rate, below which no broad lines are 
formed. The permitted interval of velocities ranges from 
« 20 000km s^^ for sub-Eddington accretion rates to 
« 1000km s^^ for Eddington accretion rates (i.e. the lines 
with FWHM > 20 000km s^^ do not exist for any mass 
of the black hole). 
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